CDT coupled to dimer matter: 
An analytical approach via tree bijections 
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Abstract. We review a recently obtained analytical solution of a restricted so-called hard dimers 
model coupled to two-dimensional CDT. The combinatorial solution is obtained via bijections of 
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causal triangulations with dimers and decorated trees. We show that the scaling limit of this model 
can also be obtained from a multi-critical point of the transfer matrix for dynamical triangulations 
of triangles and squares when one disallows for spatial topology changes to occur. 
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CDT AND MATTER 

The idea of Causal Dynamical Triangulation (CDT) [1] follows from a long tradition 
of discrete approaches to quantum gravity in which the path integral over geometries is 
| regulated by approximating it by a sum over triangulations composed of flat triangles or 

more generally polygons. The regularisation is removed at the end of a calculation by 
finding an appropriate scaling limit. In the case of Dynamical Triangulation (DT) [2] and 
CDT the triangles have fixed length sides of size equal to the cut-off and all geometry 
; is encoded in how they are glued together. In particular curvature is localised at the 

points where multiple triangles meet. This means the sum over triangulations directly 
implements a sum over discrete geometries in contrast to summing over metrics, in 
which one must mod out the group of diffeomorphisms. The difference between DT 
and CDT is in the class of triangulations summed over, with CDT introducing an extra 
causal constraint corresponding to a preferred time-slicing with respect to which no 
spatial topology change may occur. This difference between two-dimensional CDT and 
DT, which are unrestricted planar triangulations, is shown in Figure 1 (a). 

Much work has been done investigating the properties of pure (i.e. withour matter 
coupling) quantum gravity defined through CDT, in D = 2, 3, 4 dimensions (see [3] for a 
recent review): For D > 2, work has been numerical in nature and produced impressive 
results demonstrating the existence of a phase in which an extended de-Sitter space-time 
exists. There have also been interesting results related to the UV properties of the theory 
with indications that in the UV limit D = 2. In two dimensions analytic calculations are 
possible and the theory has been reformulated as a matrix model [4] showing that the 
D = 2 can be interpreted as a string worldsheet theory. The simplification for D = 2 is 
due to the topological nature of the Riemann-Hilbert action, which leads to an action 
that only depends on the space-time volume weighted by the cosmological constant A. 
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The DT and CDT programs therefore replace the gravitational path integral by, 

iF(A,T) = I dST [ ^ge-^s+SsMr) ^ $e(x,y) = TV -L e - XnT+ST ^ (1) 
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where g is the space-time metric, A g is the continuum volume of the space-time, 2? 
represents the matter fields and S g (T) the action of the matter coupling which depends 
on a set of coupling constants labelled T. On the discrete side, the integral over metrics 
is replaced by a sum over triangulations T (including a symmetry factor I/C7-) and 
the integral over continous matter configurations by a sum over discrete matter 
configurations x. Further, A g is replaced by the number of triangles nj of T and the 
coupling constants A and 7 are "bare" analogues of A and T. 

Matter coupled to CDT has been studied numerically in [5] with indications that the 
back-reaction on the geometry is not as severe as in DT This has prompted suggestions 
that conformal matter with central charge c > 1 may be able to couple to CDT [6]. Some 
analytic work has also been done [8] in which new scaling limits were found however 
this work had the disadvantage that it was unclear what conformal field theory should 
describe the matter sector in the continuum limit. 

A simple matter model coupled to CDT which can be studied analytically are so- 
called hard dimers (see [7] for dimers coupled to DT). In this model dimers with a 
fugacity ^ can join two adjacent triangles and they are hard in the sense that a single 
triangle can only contain a maximum of one dimer. The analytical derivation of the 
partition function reduces to a combinatorial problem of finding the number jV{n\ k) of 
causal triangulations (of a sphere) with n triangles with k dimers, since 
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where z = is the fugacity of a triangle while ^ is the fugacity of a dimer as introduced 
above. From the point of view of statistical mechanics, jV{n\k) can be thought of as the 
micro-canonical ensemble while iF(z, <^) is the canonical ensemble. 

In what follows we review some recent results by the authors [9] regarding a combi- 
natorial solution for determining jV{n\k) for a restricted class of hard dimers coupled 
to two-dimensional CDT. Simultaneous work by Ambj0rn et al. [10] shows that the ma- 
trix model formulation of pure two-dimensional CDT [4] can been extended to a higher 
multi-critical point in agreement with the scaling limit of the combinatorial solution, thus 
providing evidence that both these approaches yield the correct continuum limit of hard 
dimers coupled to two-dimensional CDT. Before reviewing the combinatorial solution 
obtained in [10] we present a new alternative approach to the problem using a scaling 
limit of a transfer matrix for CDT with squares reminiscent of the higher multi-critical 
point in the matrix model [10] and the related peeling procedure [11]. 

A TRANSFER MATRIX FOR CDT WITH SQUARES 

The partition function of two-dimensional pure CDT on a cylinder was first ob- 
tained using transfer matrix techniques [1]. This amounts to determining the quantity 



Figure 1. (a) Here we show an example of a two-dimensional CDT and DT. The exit boundary is the 
one marked by a dotted line. The pieces beneath each diagram represent the pieces that can be added to 
the exit boundary to construct the transfer matrix, there W represents a disc function. The final square 
piece in brackets does not appear in the pure case. Here we have coloured the triangulations by light and 
dark grey to represent alternate applications of the transfer matrix, (b) The mapping between squares and 
hard dimers. This mapping does not respect the causal structure of CDT. 

Jf{n,li,lf,i), i.e. the number of triangulations of a cylinder with n triangles, where we 
have added the boundary condition that the entrance and exit loops of length /, and If are 
separated by a distance/time t. The derivation in [1] made explicit use of the preferred 
time-slicing by identifying the edges at a constant distance from the entrance loop, or 
equivalently constant time, with the subsystem acted on by the transfer matrix. A similar 
computation in DT was performed earlier [2]; here however there is no time-slicing, 
instead the transfer matrix acts on subsystems defined as the triangles whose vertex in 
the dual graph is at constant dual-distance from the entrance loop. This difference is 
illustrated in Figure 1 (a). It is important to note that one has a composition law, 

Z(z, h,l f ,t + 1 ) = £ Z(z, /, t)Z(z, I, If, I). (3) 
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which allows one to identify the transition matrix as Z(z, k, If, 1). Technically it is easier 
to work with the generating functions ^(z,x,y,t) for these quantities. We then have the 
composition law, 

---3f{z,x,w-\t)T{z,w,y,\). (4) 
2kiw 

To compute the transfer matrix 3f(z,w,y, 1) we must construct the generating function 
for all possible ways one can obtain an exit loop from a given entrance loop in one 
step. We show in Figure 1 (a) the possible pieces of triangulation from which the DT 
transfer matrix can be built. It is important to note that in DT the effect of a change in 
the topology of a slice is contained entirely in the local moves in Figure 1 (a). 

Analogous to the construction in [4] and [12] it is suggestive to propose that forbid- 
ding spatial topology changes with respect to the spatial slices in the DT transfer matrix 
construction produces triangulations in the same universality class as CDT. Assuming 
this proposal for a moment, we may then go further and add squares to the transfer ma- 
trix as was done in [13]. The intention of [13] was to study DT coupled to hard dimer 
matter; here, by forbidding topology change we expect to recover CDT coupled to some 



form of matter. Explicitly we have, 



where the first factor adds a single marked triangle (see [13] for details) and K = zxy 2 + 
zx 2 y + z 2 £,xy i , represents the possibility of adding the first, second or fifth piece under 
the DT diagram in Figure 1 (a). Note that we have set W = 1 in each of these pieces to 
prevent topology change. Finally, by setting ^ = we can return to a pure triangulation. 
To remove the regularisation we look for a scaling limit by making the ansatz x = x c e~ aX , 

y = x~ l e~ aY and z = z c e~ a A where k is a constant and letting a —> 0. To find the critical 
point (x c ,z c ) we demand that to leading order, $?(z,x,y, 1) = cT x (X + Y)~ 1 + 0(1); this 
is equivalent to requiring jV(n, //, //,0) = 8i h i f . This requirement gives z c as a function 
of % . By substituting the scaling ansatz into (4) we obtain in general, 
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where the expression in the square brackets corresponds to 3?(A,X, £, 1), r\ is a function 
of k, F is some unspecified function and the dots represent terms which contain no poles 
in £. By closing the contour and applying the residue theorem we obtain a differential 
equation of the form, 



dr&(K,X,Y,T) = dY G F (A,Y)#(A,X,Y,T) 



(7) 



where Gp is a function derived from F and T = cfit is a continuum distance/time. We 
now present two distinct scaling limits; 

• For E, = 0; x c = 1 and z c = 1 /2. We also have r\ = 1 with = F 2 — A. This is 
precisely the result derived in [ 1 , 4] for pure CDT This justifies our earlier proposal. 

• For £ ^ 0; we find we must also set % to a critical value of £ c = — | to obtain a new 

scaling limit. We then have x c = Z c = 3 2 . We also have 17 = 2 with = 7 J — A. 
We can also extract a scaling exponent from Jlogz/J^ ~ (^ — ^ C ) CT , finding a = \. 
This is precisely the result derived in [10, 9] for multi-critical CDT. Please consult 
[11] for a discussion of the correct boundary condition to impose on (7). 

There are two important observations to make about the £ 7^ case. Firstly, we have 
only allowed squares to be added by attaching them along one edge. In the DT case [13] 
squares were added in every manner; if we were to do this here we would not find a new 
scaling limit. Secondly, the scaling dimension of z and t differs from the £ = case, 
implying that the geometry of the multicritical phase is in some sense fractal. Both of 
these properties show that the multicritical phase is more delicate in CDT than DT. To 
address why this is the case, we recall why we expect the addition of squares to lead 
to a scaling limit corresponding to the addition of hard dimers. A hard dimer model 
adds objects called hard dimers to a triangulation, with the property that they occupy 




Figure 2. (a) A bijection between hard dimers in CDT and labelled trees. The tree is shown in thickened 
lines and the labels are added using the rules on the right, (b) An example of the bijection between distinct 
trees with type-3 dimers and those without. 



two adjacent triangles and no two dimers can occupy the same triangle. The mapping 
between dimers and squares is shown in Figure 1 (b); note that in CDT the two possible 
dimers corresponding to the same square do not always respect the causal structure and 
hence it is clear this mapping breaks down. 

CDT WITH RESTRICTED HARD DIMERS 

We now review the results of [9] and argue the above model does indeed correspond to 
hard dimers on CDT. We compute JV (n, k; t), where k is the number of dimers via the 
bijection to labelled tree graphs shown in Figure 2 (a). Bijections between (C)DT and 
trees have been known for some time (see [14] for DT and [15] for CDT) so it is no 
surprise that dimers can by adding by labelling the tree. The difficulty lies in counting 
such trees since the bijection removes links and hence local interactions in the CDT can 
become non-local in the tree, thereby preventing the trees being defined recursively and 
hence enumerated. 

In [9] we strike a balance; we include enough non-local interactions to see new critical 
behaviour but not so much that the calculation is impossible. This can be seen in Figure 
2 (a); we only allow a restricted class of dimers to appear. Note that the type-3 dimers 
are non-local, however they can be accounted for by noting there exists a bijection 
between trees containing such dimers and distinct trees in which they are replace by 
type-2 dimers. This bijection is shown in Figure 2 (b). Dimers of type one and two 
only affect dimers which are adjacent to them in the tree. Hence trees containing only 
these can be counted recursively. Hence we obtain a recursive formula for the generating 
function Gd(z^) of <Ad(n,k), where ,yfd(n,k) is the number of trees with n vertices, k 
dimers and whose root vertex is labelled by d; 

G ° ( ^ )= l- Z 3(G + U+2jG 2 ) md G2fc«) = Z 2 (G„ + 2 5G2)Go . (8) 

where d = implies no label, and G\ = Gq. Note the factor of two with each G2; 
this accounts for the type-3 dimers. Solving for Go we obtain the partition function 



3f(z, £, ) — Gq for dimers coupled to CDT on the sphere. From this we can compute 
G = \. Furthermore we find in the scaling limit with % = (^critical that Go — ^critical ~ A 1 / 3 . 
This agrees with the computations in [10] and with those above. 

CONCLUSIONS 

We reviewed a recent work by the authors [9] in which a combinatorial solution of 
a restricted class of hard dimers coupled to two-dimensional CDT was obtained. The 
scaling limit of this model arguably lies in the same universality class as the full model 
(i.e. unrestricted class) of hard dimers coupled to two-dimensional CDT. This continuum 
model can also be obtained from a higher multi-critical point of the CDT matrix model as 
was shown by Ambj0rn et al. [10]. Inspired by the derivation of the higher multi-critical 
point in the CDT matrix model [10] and the corresponding peeling procedure [11], we 
proposed an analogous solution through a higher multi-critical point of the DT transfer 
matrix with spatial topology changes suppressed. This provides us with a different and 
interesting perspective on the problem. 
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